ABSTRACT
Recently, molecular covalent docking has been extensively developed to design new classes of inhibitors that form chemical bonds with their biological targets. This strategy for the design of such inhibitors, in particular boron-based inhibitors, holds great promise for the vast family of ß-lactamases produced, inter alia, by Gram-negative antibiotic-resistant bacteria. However, the description of covalent docking processes requires a quantum-mechanical approach, and so far, only a few studies of this type have been presented. This study accurately describes the covalent docking process between two model inhibitors - representing two large families of inhibitors based on boronic-acid and bicyclic boronate scaffolds, and three ß-lactamases which belong to the A, C, and D classes. Molecular fragments containing boron can be converted from a neutral, trigonal, planar state with sp2 hybridization to the anionic, tetrahedral sp3 state in a process sometimes referred to as morphing. This study applies multi-scale modeling methods, in particular, the hybrid QM/MM approach which has predictive power reaching well beyond conventional molecular modeling. Time-dependent QM/MM simulations indicated several structural changes and geometric preferences, ultimately leading to covalent docking processes. With current computing technologies, this approach is not computationally expensive, can be used in standard molecular modeling and molecular design works, and can effectively support experimental research which should allow for a detailed understanding of complex processes important to molecular medicine. In particular, it can support the rational design of covalent boron-based inhibitors for ß-lactamases as well as for many other enzyme systems of clinical relevance, including SARS-CoV-2 proteins.
ABSTRACT
MMP-9 plays a number of important physiological functions but is also responsible for many pathological processes, including cancer invasion, metastasis, and angiogenesis. It is, therefore, crucial to understand its enzymatic activity, including activation and inhibition mechanisms. This enzyme may also be partially involved in the "cytokine storm" that is characteristic of COVID-19 disease (SARS-CoV-2), as well as in the molecular mechanisms responsible for lung fibrosis. Due to the variety of processing pathways involving MMP-9 in biological systems and its uniqueness due to the O-glycosylated domain (OGD) and fibronectin-like (FBN) domain, specific interactions with its natural TIMP-1 inhibitor should be carefully studied, because they differ significantly from other homologous systems. In particular, earlier experimental studies have indicated that the newly characterised circular form of a proMMP-9 homotrimer exhibits stronger binding properties to TIMP-1 compared to its monomeric form. However, molecular structures of the complexes and the binding mechanisms remain unknown. The purpose of this study is to fill in the gaps in knowledge. Molecular modelling methods are applied to build the inhibitory and non-inhibitory MMP-9-TIMP-1 complexes, which allows for a detailed description of these structures and should allow for a better understanding of the regulatory processes in which MMP-9 is involved.